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FORTRAN PROGRAM FOR CALCULATING VISCOUS FLOW 


ON A BLADE-TO-BLADE STREAM SURFACE OF A TURBOMACHINE 


SUMMARY 

A method for analyzing the nonadiabatic viscous flow through 
turboraachine blade passages is presented. The field analysis is 
based upon the numerical integration of the full incompressible 
Navier-Stokes equations, together with the energy equation on 
the blade-to-blade surface. A Fortran IV computer program has 
been written based on this method. The numerical code used to 
solve the governing equations employs a nonorthogonal boundary 
fitted coordinate system. The flow may be axial, radial or mixed 
and there may be a change in stream channel thickness in the 
through- flow direction. 

The program input consists of the configuration of the stream 
channel annulus Sj^, the blade-to-blade stream channel weight flow, 
the inlet flow conditions, and outlet flow angle, and the 
rotational speed of the machine. The output includes the distri- 
bution of the stream function, the vorticity, the static pressures 
within the blade passages, and the variation of meridional and 
tangential velocity components from blade-to-blade and from the 
inlet of the machine to its exit are moreover generated. The 
program also has the capability to generate the temperature 
distribution within the blade passages in case the blade cooling 
is under consideration. 

This report includes a complete description of the inputs 
required for two Fortran IV programs. The first program considers 
laminar flows and the second can handle turbulent flows. Numerical 
examples have been included to illustrate the use of the program, 
and to show the results that are obtained. 
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INTRODUCTION 


An attempt is made to demonstrate the feasibility of obtaining 
viscous flow details within turbomachine passages by appropriately 
combining several blade-to-blade viscous flow solutions. Each of 
these solutions is obtained through the numerical integration of 
the full incompressible Navier-Stokes equations over a predetermined 
stream surface that extends between the blades. The set of stream 
surfaces required for the analysis are themselves generated from the 
solution of the nonviscous version of the Navier-Stokes equations. 

This report is presented in two parts. Part I includes a 
description of the input to the laminar viscous flow Fortran IV 
program and the interpretation of its output. A numerical 
example has been included to demonstrate the use of the program. 

The example includes flow in a radial inflow turbine. 

Part II presents the procedures for computing turbulent 
viscous flow details within turbomachine passages. This part 
discusses the input and output, with a numerical example 
to demonstrate its use as a turbulent viscous flow analyzer. 

The present report forms the final volume of a three volume 
report. This volume is organized as a user's manual to serve an 
individual desiring to use the program. For those who wish to 
study the calculation procedures and the mathematical analysis 
in detail and need to modify the program to suit their needs, 
this report will have to be used in conjunction with the theore- 
tical procedure which forms the contents of Vol. I and Vol. II [2, 
7]. The entire report is organized such that a design engineer 
desiring the use the program need only read the sections 
NUMERICAL EXAMPLE and DESCRIPTION OF INPUT AND OUTPUT, of Part I 
alone or both Part I and Part II. If only laminar viscous flow 
details are needed then reading Part I would suffice. However, 
if turbulent viscous flow details are sought then both Part I 
and Part II should be read. Both sections have been organized 
into decimally numbered subsections, and accordingly a decimally 
numbered system is used for the equations and figures. This 
numbering system is adopted for the ease of handling the two parts 
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with overlapping information and frequent cross references. In 
this decimal numbering system# the first number refers to the 
part of the report in which the section# equation or figure 
exists and the subsequent numbers refer to subsection. 


PART I 


1.1 NUMERICAL EXAMPLES 

To illustrate the use of the program and to show the type 
of results which can be obtained# a numerical example is given. 

Radial Inflow Turbine ' 

An example turbine rotor profile is shown in Fig. 1.1. First 
a meridional plane analysis was made by the quasi-orthogonal 
method of reference [1] to establish the stream channel geometry 
and the stream annulus thickness distribution. The flow patterns 
are investigated on this blade-to-blade stream channel# S^^# 
located midway in the passage depth of the rotor as shown in 
Fig. 1.1. The results presented here are case 2 of reference [2]. 
The primary reason for the selection of this specific rotor is 
that a substantial amount of experimental data is available for 
comparison [3]. The input data for the program are given in 
Table l.I. A sample of the output of the program is presented 
in Tabale l.II. The streamlines are designated by the stream 
function ratio ’1^/^ total’ streamline contours are plotted 

in Fig. 1.2 along with the streamlines determined experimentally. 
Also presented here are the velocity profile at three radial 
locations as shown in Fig. 1.3. For the purpose of comparison 
the inviscid flow profiles# indicated by dotted lines, are 
superimposed on the viscous profiles in Fig. 1.3. For a more 
detailed discussion of these results# Vol. I of this report should 
be referred [2]. 


Radial Compressor Rocor 

This example is a rotor of a radial bladed compressor of 
reference [4]. The rotor profile is shown in Fig. 1.4. First, 


3 


a mid-channel stream surface, S^, and the stream annulus thickness 
are calculated by the meridional flow analysis of reference [1] . 

A typical distribution of the flow properties on the blade-to-blade 
stream annulus, S^, of Fig« 1.4 are calculated and the corresponding 
results are presented. 

The results of radial compressor rotor flow are presented 
solely for the sake of demonstrating the capability of the 
program. Therefore only the final results are discussed. The 
preparation of input and interpretation of output is similar to 
the previous turbine example. A comparison between the predicted 
streamline contours and those determined experimentally in 
reference [4] are plotted in Fig. 1.5. The experimental evidence 
was obtained by tracing photographs of streak lines from the 
rotor segments under the same operating conditions reported here. 

A good agreement is generally observed. The predicted meridional 
velocity profiles across the rotor passage at three different 
radial locations are • illustrated in Fig. 1.6. Shown also in the 
same figure, are the calculated velocity distribution using the 
inviscid blade-to-blade analysis of reference [5]. 

1.2 DESCRIPTION OF INPUT AND OUTPUT FOR LAtllNAR VISCOUS FLOW PRO GRAM 

The computer program requires as input appropriate gas con- 
stants, operating conditions such as inlet temperature and 
density, inlet and outlet flow angles, weight flow, rotational 
speed, and information regarding the overlaying of grid in the 
flow passage. Figure 1.7 shows the m-9 coordinate system for a 
typical blade-to-blade surface which includes the physical flow 
domain. Also required as input is the stream surface geometry 
and transformation parameters at every grid point. The latter 
part of the input is obtained as a consequence of applying 
Thompson transformation to the blade. Ovitput obtained from the 
program includes velocity magnitud and direction at all mesh 
points (meridional and tangential velocity) , value of nondimens ional 
vorticity, Reynolds number and pressure coefficients at each mesh 
point. 


Input 


Figure 1.8 shows the input form indicating the variables as 
they are punched on the data cards. There are two types of 
variables > geometric and nongeometric. The geometric input 
variables are shown in Figs. 1.9 and 1.10. A description of 
all the input variables is listed in this section. 

Further explanation of key variables is given' in the section, 

I 

Instructions for Preparing Input. The input variables, in order 
of their appearance, are as follows; 


AMU 

PRDNO 

TIP 

TBLADE 

RHOIP 

WTFL 

OMEGA 

BTI 

BTO 

BLADE 

BIN 

RTIP 

RIN 

MXBI 

MXBO 

MX 

NBBI 

ITRB 

IPRINT 


2 . 

Kinematic viscosity (N-s/m ) 

Prandtl number 

Total inlet temperature of the turbine (*K) 

Blade temperature (°K) 

Total inlet density (Kg/m^) 

Mass flow rate per blade,, per stream channel, (Kg/sec) 
Rotational speed of the machine (revolutions/sec) 

Inlet flow angle to the machine (deg), (see Fig. 1.9) 
Outlet flow angle to the machine (deg), (see Fig. 1.9) 
Number of blades 

Inlet stream channel thickness, (m) 

Tip radius of the rotor, (m) 

Inlet radius of the machine at upstream boundary AN, (m) 
Number of vertical mesh lines from AN to BM inclusive, 
(fig. 1.10) 

Number of vertical mesh lines from AN to FI inclusive, 
(Fig. 1.10) 

Total number of vertical mesh lines in m-direction 
from AN to GH; maximum of 50, (Fig. 1.10) 

Number of mesh spaces on cj) direction between AB and NM; 
maximum of 30 (Fig. 1.10) 

-1 for turbine case; +1 for compressor 

The number of iterations after which the output 

should be printed. 
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0 


ISPL 


AMM(I) 

R(I) 

SAL (I) 
BRO(I) 

UTR(I,N) 

VTR(I,N) 




UIC(I,N) 
THETA ( I, N) 


For a straight blade, ISPL is assigned any integer 
value other than zero, otherwise it is left blank or 
zero. If ISPL is zero, then appropriate number of 
transformation parameters follow the stream surface 
configuration 

Array of m coordinate corresponding to mesh points (m) 
Array of r coordinate corresponding to AMM(I) array, (m) 
Array of sina, (see Fig. 1.10) 

Array of stream channel normal thickness corresponding 
to the AMM(I) array, (m) 

Transformation parameter 

Transformation parameter 

Transformation parameter 

Transformation parameter 

Array of <{) coordinate at 


(= 1^) [Note: 3X = 

(= -^) [Note: dY = d(J)] 
(= 

{= 

each mesh point. 


Instructions for Preparing Input 

Units of Measurements : 

The International System of Units is used throughout this 
report. The program does not use any constants which depend on 
the system of units being used. Therefore any consistent set of 
units may be used in preparing the input for the program. 

However, the angles must be supplied in degrees ;inless indicated 
otherwise, and the rotational speed is given in revolutions per 
second. Since any consistent set of units can be used, the 
output is not labeled with any units. In the examples supplied, 
the International System of Units are used. 

Stream Surface Geometry : 

The stream surface on which the flow is being determined, 
is a surface of revolution. Hence it is sufficient to define the 
mean stream surface of revolution (as seen in meridional plane, 
Fig. 1.9a). The stream surface required for the analysis is 
generated from the solution of a nonviscous version of the Navier 
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Stokes equation as suggested by Wu [6] . The mean stream surface 
is defined by a set of m--coordinates and corresponding arrays 
of radii/ sine of the angle made by the m-coordinate to the 
machine axis and the stream thickness for each 'm* . 


Inlet and Outlet Fl ow An gles : 

The value of given as an average value on AN. The 

exit flow angle is found iteratively. Normally one or two 

trial runs of the program are required to arrive at the exit flow 
angle. The iterative procedure is as follows. Estimated exit 
flow angle, ^@xit' used to specify the values of 

stream function derivatives in the m-direction through the 
following relation 


^ _ _1 

3m “ 2 it/z 


tanB 


exit 


exit 


( 1 . 1 ) 


The flow field equations are then solved for the boundary function 
ijj given by equation (1.1) to Obtain the velocity and pressure 
distribution throughout the stream annulus of Fig. 1.7. An 
evaluation of the torque developed by the annulus is obtained 
through the integration of the pressure and shear forces acting 
on the blade surfaces. The change in the annulus momentum between 
the known inlet and the estimated exit flow conditions is deter- 
mined. If now the value of the predicted torque is not equal 
to the rate of change of the angular momentum, then the direction 
of the exit flow velocity is altered. The whole procedure is 
repeated until a satisfactory result is obtained. 

Defining the Mesh t 

A finite-difference mesh is used for the solution of the 
basic differential equation. A typical mesh pattern is shown 
in Fig. 1.11. The mesh spacing and the extent of the upstream 
and downstream regions are determined by the values of MXBI, 

MXBO and MX of the input. The program is presently set up to 
accept a maximum of 30 for MX. This is not a necessary constraint, 
the limit can be changed with very little effort. The values 
of MXBI, MXBO and MX should be determined such that the mesh 
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which results has blocks which are approximately square. To 
achieve this, a value for NBBI is first chosen arbitrarily 
(15 to 20 typical) . NBBI is a number of mesh spaces spanning 
the blade pitch S, where S = 2 tt /BLADE. Dividing S by NBBI 
gives the mesh spacing DN in the 0-direction in radians. Multi- 
plying DN by an average radius (RIN) of the stream channel gives 
an average value of the actual mesh spacing in 9-direction. 

The value of axial chord of the blade should then be used with 
this tangential mesh spacing to calculate the approximate 
number of mesh spaces along the blade in the m-direction. This 
will give MXBO once MXBI is chosen. 

Format for Input Data ; 

All the numbers on the first and second cards (Fig. 1.8) 
beginning with AMU and BTI respectively are real nxambers (punch 
decimal point) in a 10-column field. The numbers on the third 
card beginning with MXBI are integers (no decimal point) in a 
five-column field (Fig. 1.8). And finally, all other numbers 
following the fourth card are real numbers specified in E format 
in a 14-column field. 

Output 

A sample output is shown in Table 1.1 for the radial inflow 
turbine rotor example. Since the complete output would be 
lengthy, only the first few lines of each section of output are 
reproduced here. Each section of the sample output in Table 1.1 
has been numbered to correspond to the following description: 

(1) The first output is an echo of the input data from 
Fig. 1.8. This comprises of sections (la) and (Ic) . Though (lb) 
is not strictly an echo of the input data, it is useful for the 
future correction of the input data as explained in the section. 
Defining the Mesh. This is a printout of DS and DN, the grid 
spacing in the m-direction and ©-direction, respectively. 

(2) This is the output which ’consists of blade transformation 
parameters. The output is displayed for each mesh point and 
begins from the line of constant n = 1 till the last n- 
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(3) These are the mean relative velocity at the inlet (WMIN) 
and the inlet Reynolds number. 

(4) This is the maximum iteration required, the maximum 
values of the stream function, and the maximum value of the 
vorticity in the flow field. 

(5) The values of the stream function, the meridional 
velocity, the tangential velocity, the vorticity and the pressure 
coefficient are displayed for each mesh point. These quantities 
are presented in groups along constant n lines. 

1.3. PROGRAM PROCEDURE 

The program is segmented into seven main parts, the sub- 
routines VORT, STREAM, STAPRD, STORE, TEMP, CONVERG and sub- 
routine function SPLINE are called by the main program. All 
the subroutines and their relations are shown in figure below. 

All information which must be transmitted between the seven 
main subroutines is placed in COMMON. 



Fig. 1.12. Calling Relations of Subroutines. 

The program was run on AMDAHL 470/V6-11 central processing unit 
with six magabytes real data storage. The program is presently 
set up to handle 30 by 40 mesh points; this was chosen to keep 
the running time reasonable. If, however, the user needs to use 
more mesh points, or if there is a storage problem on the user's 
computer, the maximtim number of mesh points can be increased or 
reduced as the case may be by appropriately altering the dimension 
and common statements. 
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PART II 


2.1. NUI^RICAL EXAMPLE 

The equations formulated in Volume II [7] are programmed 
for numerical solution using the finite difference technique. 

The program is arranged to handle general flow within turbo- 
machinery/ which may be of the axial, radial or mixed flow 
type. To illustrate the use of the method of analysis a sample 
problem will be discussed. The input and output of the program 
will be presented. 

Nonadiabatic Flow in a Mixed Flow Turbine 

The versatility of the present method of solution is pre- 
dicting the behavior of the viscous nonadiabatic flow in turbo- 
machinery is demonstrated by analyzing the flow in a cooled 
mixed-flow turbine rotor. The rotor geometry is shown in 
Fig. 2.1 together with the shape and thickness distribution 
of the mid-channel stream annulus, The blade-to-blade shape 

in the physical domain as well as the boundary fitted coordinates 
employed in the solution are shown in Figs. 2.2 and 2.3. 

Additional summary data related to the blades geometry and the 
configuration of stream annulus can be found in Table 2.1. 

The following operating conditions which correspond to the design 
point of the turbine are used in the analysis. The program out- 
put is presented in Tables 2.2 and 2.3. 

Turbine inlet total temperature, T. , 1083®K 

^ 3 

Turbine inlet density, p, 1.0060 Kg/m 

Mass flow per blade for the stream annulus, M, 0.00320 Kg/sec 

Rotational speed, f2, 38,500 rpm 

Inlet flow angke, (see Fig. 1.10) , 62.5° 

Exit flow angle, (see Fig. 1.10), 63° 

Meridional component of the relative velocity at rotor inlet, W^, 
66.2 m/s 

Flow Reynolds number at rotor inlet. Re (= — — --P ) , 3.34x10^ 

Prandtl number, 0.8 inlet 

Working fluid, air. 
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The results are presented in Figs. 2.4 through 2.9 as 
contour plots for the distribution of the stream function, the 
velocity, the kinetic energy of turbulence, the static pressure 
as well as the temperature within the blade passages. The 
detailed discussion of these results will be found in 
reference [7]. 

2.2. DESCRIPTION OF INPUT AND OUTPUT FOR TURBULENT 

VISCOUS FLOW PROBLEM 

The program requires as an input the configuration of the 
stream channel annulus S^^, the inlet flow conditions, the rota^ 
tional speed of the machine, and the blade geometry. Recalling 
that all the flow calculations are carried out in the unit 
square of the transformed domain, therefore the blade input 
geometry is supplied to the program in the form of the trans- 
formation parameters 5, B, Y and J. As pointed out in 
reference [1] , these parameters may be specified for any blade 
geometry using Thompson code for the automatic numerical 
generation of boundary fitted coordinate system [13] . 

The program output consists of the distribution of the 
stream function, the vorticity, the kinetic energy of turbulence, 
its dissipation rate, and the static pressures within the blade 
passages. The variation of meridional and tangential velocity 
components from blade-to-blade and from the inlet of the machine 
to its exit are also generated. In cases where blade cooling 
is considered, the program has the capability to generate the 
temperature distribution within the blade passages. In order 
to keep the computer time within reasonable limits, (usually 
less than 5 minutes on an ANDAHL 470) , the flow domain has been 
divided, for all calculations, into 30 step sizes in n direction 
and 40 in the i; direction, with the greater number of nodes 
distributed in the meridional direction. 
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Input 


The description of the input variable, the input form • 

(Fig. 1.8) and their formats are presented in Part I of this 
report and are valid for the turbulent program. Also the 
comments made about Units of Measurements, Stream Surface 
Geometry, Inlet and Outlet Flow Angles, and Defining the Mesh, 
in the section, "Instructions for Preparing Input," of Part I 
hold good for the turbulent program presented here in Part II. 

Output 

However, the output for the turbulent program needs some 
more explaining. The saunple output is shown in Table 2.2 and 
Table 2.3. Here again, since the complete output would be 
lengthy only parts of all relevant output are reproduced in the 
two Tables 2.2 and 2.3. Each section of the sample output 
in Table 2.2 and Table 2.3 are sequentially numbered in sub- 
divisions to correspond to the following description; 

(1) The first output is an echo of the input data. This 
comprises of sections (la) and (Ic) . The output in section (la) 
is identically described in Fig. 1.8. Section (Ic) displays 
the blade geometry data. Though section (lb) is not an echo 

of the data, it is an important derived quantity DS and DN, 
the mesh spacings. These are displayed to help adjust the mesh 
size. 

(2) This is the output which echoes the blade transformation 
parameters. The output- is displayed for each mesh point, and 
begins from the line of constant n = 1 till the last n. 

(3) These are the mean relative velocity at the inlet, 
the Reynolds number at the inlet, and two intermediate para- 
meters . 

(4) This is the maximum iteration required, the maximum 
value of vorticity, the maximum value of the stream function, 
the maximum value of kinetic energy of turbulence, and the 
maximum value of dissipation function. 

(5) Finally, the values of stream function, meridional 
velocity, tangential velocity, vorticity, kinetic energy of 
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turbulence, dissipation energy, Reynolds number and pressure 
coefficient are displayed for each mesh point. These quantities 
are presented in groups along constant n lines. 

2.3. PROGRAM PROCEDURE 

The progreun is segmented into ten main parts, the s\ab- 
routines, VORT, STREAM, ENERGY, DSPAT, TEMP, STAPRS, STORE, 
CONVRG and subroutine function SPLINE are called by the main 
program. All the subroutines and their relation are shown in 
figure below. All information which must be transmitted between 
the ten main subroutines is placed in COMMON. 



Fig. 2.10. Calling Relation of Subroutines 




TAsLt 2.1 BLkDZ GEOMETRY DATA FOR T.HE 
MIXED FLOW TURBINS ’ ; 


Blade Upper Surface Data 


Meridional 
Distance m(m) 

0. 02S39J- 03 
0* doCOOU- C2 
0. 1 eOCOO-Cl 
0*4!35CCO-C1 
0.200aOU~C 1 
O.S«2dOO-Cl 
0*3 904 00'- Cl 
0. *62300-01 
0.a3Vft0O-01 
0.U1970J-01 
0.6701 10-01 


Angular 

Coordinate 

<> 

0.5C273 
0 .SovOO 
0.&6a6C 
0 .56220 
0.56J5C 
.0 . 5*^00 
0 . s<»e •« 0 
0 .515*0 
0 . * 2 o 0 0 
0 .23660 
0. 23a^bO-01 


Derivative 


'-0. *5676 
— 0. *6035 
-0. *ei 20 
-0. 35261 
-0.3*501 
-0. 520 1*0-01 
-1 .7545 
-7.61 16 
-15.V62 
-33.6** 
-30.7*3 


2nd 

Derivative 


-7.255* 
6.359* 
-1 1.995 
*6.208 
-*3.5*2 
1 5*. o6 
-813.31 
-926.77 > 
-12*2.5 
-3206.* 
-3 *v7.6 


Blade Lower Surface Data 


Meridional 

Angular 

Coordinate 

Derivative 

2nd 

Derivative 

Dist 2 ince Ri(m) 




0. * 2S3v.>- 03 

0 .6*7100- 02 

0.*S673 

5.7000 

0.6O0C0O-C2 

0. 122000-01 

0. *6655 

-3.2*67 

O. 1 60 00 0— C 1 

0 .160000-01 

0.46729 

0 • 7*668 

0.336000-01 

0 . 1 O900D-01 

0.40262 

- 1 5.325 

0.290300-Cl 

0 .20<}C0O-01 

0.322 13 

-13.767 

0 . 3 *t 2 6 Oo-" C 1 

0 .21 5COU-0 1 

-0. 26520 

-20V.96 

0.39S*00-Cl 

0. loOCOU-Ol 

-2 . 06*3 

-* fO. 30 

0.*623 0i>-01 

-0. 116000-01 

-U.6949 

-917.03 

0.639*00-01 

-0 ."ia 2 0 00- 01 

- 16.993 

-1*95.0 

0.6 J57 0D-C 1 

-0.2791C 

-30.3*2 

-2076.9 

0 . c 6 3 63 0- 0 1 

-C.530 7C 

-60.7*3 

-*303.7 
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TABLE 2.2. INPUT DATA CARD LISTING 


e> 

c 

^ c 

c t 

u .‘0 
s m 
u o 


c zc 

c:i c 
_/C o 
u. « o 
O 

» t\> 


i 

c 

!r. 

D 


e 

C.OH c 
— c a* o 

cc 

in 

z ir^ 

< 

a? 

c 



p 

i 

• 

i 

r 

1 

c 

i 

1 

c 

o 

L'O Z C 

CO- 

>o 

< c z 


• 

' N 

so 

c 

KC 

c 


1 

’ U 

c a c 
o 

z c -p c 
-c- *r o 

H o 

■/, c • 

— f \‘ — 

j ft — 

c 

c 

< 

u 


S J £c 

< i O Ci. Vf . 


C £. 


z o 

sca- 



•C 1 


it • 

r.w- . 


o ^ 

«- i 

o 

i 

1 

<- 

t 

M 


I- 

C- »*' 

M 

i 

— ■ 

zo 

4 

z 

c • 

*5 


o 

1 

X c^ 


n. 

sr 


<? 

o 


o 

o 

s 

• o 

o 

in 

Z 3 

— o r>c 

,m 

r 0 - H- c a r. {/> n 

.< N 

a «x 

C V 


<f z 

OJ 

’ f 

in 

o 


1 - -y 

* 

; i 
' 1 

X • 

L. , 

o 


nr 

fd 


rH 




H 


Nf^\OocarriOcccO'd^CotrN^^cr5<fotf)^.‘^vOr^inc'd‘6ccinocJom^sC® 
O O vO UO C so O' <\Ht LO vC IT wD N in CO ^ 'C* CD p-» C* CV ^ -3“ sC CVJ vC N C 00 

Qia‘*0'C^O0'Ccoo*^*^«^^-i-*^^.^(VfCj0Jr0ini^r0<t<-^<'ininin*nv0v0v0vCs0N0s0 

COOOOOp - ‘ ^ ^ 



J C O O < CM o CJ N rj m 'C* — CO N c <3* O O OvJ X' a' ^ O CM C in 'O O IT O' <t lO a* IH IT; CJ fT:- 'C 

< X cr o; X fO'X c m N o sC r/ N c c m h- N o N c: 0^ o N a: O’- r> CD X O' r CM O' vo c 

X CM Csi ^ 0 ‘ N N o O' N ^ c O' ^ O' m ro ft'; in CM X o* ^ 3D c o N o in O'-. r; in N ••t 3 j in 

<?'«-of^in-<TCN^oooo^ocxc CO c <r w o c o c o X in c o in X n c 
opooxNxr/Ovo^'d'®r*'^f^W'C^J*'^^*^o-XNNOco-^oicxr)Xin'TncM-- 
C C O O O O O O O © X N O X in -5f fO <VJ OJ OJ ^ ^ CM C'-* CM cv; CM ^ ^ c c o o c o 

♦ V •' » •♦ ? • • ^ • • • • » • * ir « • ••••»(►• • • • « *••••*• • 
- . o r. c c o c. o c c c 6 c o c c o c c 6 o o o c o e c c c o c c o c^ c o o 


V) 


I i t 


I I f 


I i * 


( t I 


I < I 


I I I I I I I 1 f I I t I M 1 


< s d. o X Ov’ o r c xc ^ c'.i c o c N x lo in n r> r> cv cw o o c oc o^o o 
o X X X N N \C“ o c c in IT in in ^ ^ <c -e <r ^ ^ <7 ^ >T <• 'sr f^' ro n n n 

C c O O O O o O o C c o c p O O c G o O C C ^ 



SOOOOOGOOOCOOCOOOOCpOCGOOCOO O.G QOOOOOOOOO 


< I 


.. ■ ' N ' ■ 




-V” 

o 


18 


OPJGiNAL PAGE IS 
OF POOR QUALTIY 


TJ 

0 ) 

:3 

c 

•H 

4 J 

c 

o 

u 


CN 

CN 


w 

CQ 

< 


l! 


1 

C. I 


\L 

<. 


u 

> 




* i 

C ^ : 

• <r 

j 

<r 

1 

j 

I c c < 

r 

i 

‘ *■ 


i c o c ■ 


, ^ 

« M M ^ 

M 

< ^ 


a. 

! OOC ; 

c 

c> 

, c c o 

C ! 

C o 1 

1 * , 

*» 

■ CC O 

c 

c 

j or o 

c 1 

! CO 1 

j . : 


C C 

c 

c 

‘ c c o 

o 

r er ' 



c c c c c c 

o c c c c c-c 

C C C c 

r c c c- c 

CCOOriCCCCC.CCCCCC 



wm 

•• •— 


«► 

O 

r. C 

I C C C c c 

«» 

o 

c r 

' oc c c c 

c 

c 

c c 

' c c c- c c 

c c 


c c 

: e r: c ^ c* 

c. r 


^ S W 1*' W W W W ^ 'V 

c * occ;er:c^c* cc 

c c C- C c o c c c c* o c c c c c c- c o C c c c o c c c* « 

COCOCCCCCCCC'OCOCCCC®C*CCCOCC‘C:I 



OF POOFJ Q<-*' 


19 


m 

CN 

03 


ro 


m 

cn 

in 


iN 

t-3 

< 

cn 


‘in 


u. 

c 

u 

• 

IT. 

u. 




CV| 

VO 

r^ 

00 

9 

lO 

VO 

II 

s 

H 

s: 

:s 


c 

z 

u. 

& 





• 




•Y 

m 


on 

5 

ro 


H 

h 


vO 

O 


4 

00 

<— ♦ 

• 

<Ji 

f-f 

O 

. V 

1 

o 

o 



o 

« 


It 

4 

o 



o 





II 


♦4 

It 



< 


2 


C/3 

2 

04 



H 

cn 



CO 

a 

t 


o 

a 

N 

r» 


o 

a 

c 


I! 

2 2 

u o 
u a 


I 


♦ 


; I i ; t : 

^MncecciciT. - -<cc•^*o-^' — cooo«>Nc«KrJir 
^ r- N c o s r cv r r. N i c. r — ? fv< «? < r- — Cw’ c n.c •** r. 

— ff<*r*. r N.r«* C*rr*<v so — r^^f»-er Kr^K^^C‘t^. n 

rj f . c o f*“ 2 * r. r <«• ss <v 1 “. N o r-j r/ m ^ c c o o o c tf. r% — c 
<vcMr^if'incNKNCCo:o-?‘C'0 0'C'0‘c<?c»cc*o*c.cc? 


o c c o c* c c 

M I ( t i I ) t * t 


1 


I 


I I I 

I 


• — — cooc o 
I t I t I t < ( I I 


J 


ryL'i^rc**-nr.ccri?^C‘e^?^o — cOk^rjNOc^rew 
r<*ir. cr'C IT. «— aocCM>iN^ tr, Nr, CCS w< 

c c fv c c c a c c -e c :w c r- c N r — N o c ^ c V* — N <• c r, 

r. evj r c c c rg c N c N K r. — c c «? N W — e O' tT/ c tniT <v N 

r<Ncaosf^^c<sc — — c—*nrsr^o*C'eroCN 

c c -e c c c O' « c rg s s ir c N - c N c c r «« r? — N r. c c < 

o N n o tfj o c- <N. c r N c o N i'> re s c -e s ir. r< O'® o N s 

N-csN'g^u“x^^ — ojre<‘Scs-ceN«'wONun^w — 

^ a C S N w r - c S' a s C c JT. c c r, r, r CN. r. ft* ftj — 


i 


i 


s* ^ r a c r. N < N p. r. cv T4 f p. <r r; c u': N — a a a- N N r- 
<• s •'*. s t/. »? n <• o D ^ o S' c- s N N \r. s r ® x r.</> n 
c s < a r o c a r s o r? u*» a* s p.’ < e c ^ r n r* r s c n s a 

a n r* fw « N ^ r o N c © N a pj r. c u*^ fu (\» c o — in r» o r*‘ r ^ 

o © c c c o — — r < s p- 'J' — c r < N o •• c (\ p- fj UV e r. ce N 

cocaoc-coccc — — Nn^cC'U*. r. <©NCCcao?u 

*-wn^cS‘'ffpirNo« 

i . ' . i ' -<vir#«u: 

' • !■ i* ' 

I ■ '' : 

. i . • » 

' ! 

a c c p- N o © cir T; K ir; a c c a S' s uo c IT N O' fg in IT- c « r 
a c '£ N s- — s e XJ a c- c N ^ \n s c « N ^{ r» c ^ p- s c 

< r * a r; pu o ^ — c r c — 4 c a r s © -e <\* c K ^ c r in N c < 

pg N a c — a c © c Pi c* <■ a N r, c r.' in Pw* a a a c N c N p^ N 

— c < c N C’ r a ui in N - N o e p. a p.* o c < c a e c p; a a 


- o c c o a 


- — p; r c a N s re IT s c IT. c c c in c © r. s N a 

— — — pjr«oN©p«in — N^c 

— -,PJPJtT>^ 


! 


I 


ON in IT 

S *- S' c 

PJ p; C IT 

< N o < 

C'C r o 


<r c 

C N 


pj ©— a r © eg a a s N a c c *n © c c « c a IT 

c P4 in a c r. s o c P-- c a — c — c o o pj c © p^ 

me c — c c c a c- © N N N c •- tn ^ — c 

<r a 'n pj in in IT a pv c o a pj N ^ © p, p.i c r. a a r 

r c N r c vM a LO — — © u* Ln c N c r> N <• e c c c g 


C N ^ 


— p> <- N — m N r, ^ — 0 p. r a c c N pj N c N c « pg c 

— — ^ — cwc^rjr. «-<'t;icc©cop.<NoinpwoocN 

— — cj Pi r>< < r. 







00 




, 

i—i 





1^ 




1 

iH 



• 

m. UP c pi p: vn p* 1 

iH 



Mt 

- c m. r ^ in a 

on 



Ui 

3 a a c s c 'n 1 

00 


0^ 

> 

©cm—©-- 

CNJ 

o 

vO 

• 

P. c' - — - 


fN 

01 

m 

- - C C C' c o 
r 

It 


H 

< 

m 

O 

• K 

Q 

o 

O 


1 ; 

♦4 

o 

• 


; j 

O 

9 

O 


; ' 1 

Z 

o 

1 


1 ; 

>1 

1 


I 

’ 5 


© cm c r. e N . 


f 

- c N a 0 1 


N \\ c © c . 

> 

c rv - c p*^ i 


a «^ © N m 1 

S 

a C N N C : 

a* 

• • • ♦ • 

£ 

•T © c r, r 1 

1 

— cv 

1 1 

f 

< 1 ! 


H ' ; 


1 


! — r.c^ 

I I I. - 
I 


ap.ma© — rm-m-r,pjoN©-^coo<-N — 

— r. p. < c pg r. c X s p. c a p- <r « c c s o m © tn 

a r- p; ^ o s a m r c — in c c c c c r c a p* c 'C N 

c © r. m a c 3 m © s r 'n c c m a r, a © a — ^ c* a 

a — a < ^ r a m © m c r* p^ ^ c — c < — ’-np.* m p. n 

r. c c m c p» a © © c — p? m N a T 'S p^n ©•- c p 

I — p: r: IT 


I I I I 


o 

50 

• 

U 

u. 

r 1 
N ‘ 

in ! 
Z 1 

LO 


t 

o : 

rH 

11 

r— 

c o ©• 

11 

&H 

m 

• • • 

i C Ol.' 

Z 


c 1 

05 

D 


t' ’ 

W 

O 



O 

1 

^ ' 

H 

H 




•ooooococoocco 


; 1 

i i ■ 

p- c c N m e r © *'- m c> m K © — N c c r P» n r c © pj r> N p; 

c c c N c a rt c — e. ►. a m © c o •. c - a o a c — m o p; c 
— c r — p‘inp'P‘N CO ar m.mr N r e m - rniK ©m© r. m c 
c©acrsmocap*cN©-rcinCNcaccaac'r-. 

(vp-rcinmcN^-Nrarao ss ©©©s©©s©©©c^ 


c o C C C 
I I I I I I 


» 


f I I t I I 


I 


♦ 


*♦11111 


f I I ♦ I I t t 


t 


a©— mp.'asscoointnmro — p.©©©o — m — ccnn 

cc in — c rmc a cr s - © c c N - N mpgcc s p- in m 

Nr.NcNCPJtnecNNc^-rNmcrr.r, c — — NC©co 
© — m c r c c o c c rg a s — — s — pg m in N N pj m r m c 
N c c o c p; © o — © m o m '-p © < in — c w r a r p n o s - c 

p. n n r r c pi r. - © m P’ a © m pj a c N m c K in c © N s c r. 
©ciniPp%-piine cm. c a p, s %c —c pj© n c m. in c c c n 

in a m pj & a c c r. — o © c a N s. c c m in c c c c c c c c c 

mpipjpjp; — — — — — — 

i' ; i ! 




f 


.c — a N r c c a c a r c m. - c r c c a c c a o N a. m c N 
r o o K m* r p. c o c c - m. c po © Pi. a a p. s a m. c m p. r, 

a c a c a © Pv‘ •" c X ^ ^ pg c — — — n s* c c c © — a c ^ — 

- a m — — s © p: © m c m N a c — r c c fj in r* a- r r c P. c ? 

c 'C N N © p; a in <■ p- r; n s a a m p > c c- p* r •- a o Pi < r s: <. 

o — occ — -pgrecce^Nr. — p. cc^Ncr,p.©p. N© 

---.pgrcmN©*“cCc--^CN 

• ! . \ fW' r,c m. m 


1 : i 

a— m. r cc©c 

a — 0 m N © pg e p © pg !« 

CPi©N0CmC--©CN 

P^PJC— P.CPJ — — NeC-, — 

^ r a p- c N c r. ^ f'** c - m c c S'N c s r Pi in p k c o 

— in c — r. c — c pj c © a 0 - — m r c c a N m p», p,- © Pi p; c r-. 

— — — — — p.' pg r c c m c a- © — r c c fn n m. © c r, o c 

— — — ^ pg p; n r r < m. c c 


r - < c a r P-* m. p m rv n © a c a © < 
•>CNm;a*"©Nrr as — N.re — r 
N p r U-. N c r c c m N a P’ c p- © 
c m, a m. m, ■ a c c c m pj © a n c m 

<• e e s. Pc ir. p. C r. */: 


I 


-Nccoo©c—mo- 
r N in N - c © m s o - r* 
r a N c a - c m >£ c r m, 

N©cN©©mCNcas 
Nc-ccp*©— p. p.^ac 


I 


s c n c a N Pi o c o p c © p- in r • 

m r* — c c a o — r c — c o p- pg ; 

c a r o © © c - c r un c p* . — p. o f 

p p p r, c c p* Pi o pj a a © N a N • 

s c r. m a — a © c 'n r* a m, n — p, i 



c n a < p*i o pj c r ©• N c r r m r n c c r. m o* p - < r. 

c* e — r a — © a m ^ — m © r cg c m a ©• p p c m •© a c c < 

n — c N c Pv c © c^. <■ a c s pv o — *c c- © a c m > c m r m? - n 

rv K c p p. c p; s ^ © r N 0 pg — c e ^ a c c o <• r n c a a m 

m r © r - c c c o c - - a c m, r- . a r c m. c r c* © c 


♦ I t I I ♦ 


o o a r-4 r m. * 


‘ © a m - m. a 
I -r 
I * 


cc 

7 


p; p.! © Pi c a c m- © — t* a c a a © c © c a a p. c c o m 

c c - c p‘ in c m © pj c c p. c vo c c PI <* © - a N o m m ©cm 

p*/ p' N N C — c: o m © c r. 0 m © c m c © ^ n s r. c ^ s r. r, © 

N r. c — o © o N a a a N c c c. ^ © — c c c a N N a c m c 

m. m © c < r j © ^ c m c c c m ©■ c c ^ c © c* r, n c o © © 






• ♦ » 


a c c c < c m m c N a © Q o -» p j p* m c c m p; © p* c m* © "n © 

- - Pi c c 
I I ’ I 

i • 1 I ♦ i < 

< 

! 

N - c c a Pi N © a PV « a a r*. © n c m in — r. c — c c c c m ' 

m c r. c m c a ^ © Mn p* © -^ p< m a -n ^ o c a © n <r c m 

2 r> f 4 ^ ^ r, f.j p: r. fp r. ^ m. m. c m m r.-- r. .■? © c © 

ooccococccoccscocoooccoccoco-'- 

UOOCOCOOCCOOCCOOCCOOCCCOCOCOOC 

U.OOOCC ecjooooaC'Oci 
C. t I 1 I 1 I i I I I i I i I t I I ♦ I * < I t I 


= I 


20 








REFERENCES 


1. Katsanis, T. , "Use of Arbitrary Quasi.-Orthogonals ' for Cal- 
culating Flow Distribution in the Meridional Plane of a 
Turbomachine/' NASA TND-2546, December 1964. 

2. Khalil, I. and Tabakoff, W. , "A Calculation Procedure for 
Viscous Flow in Turbomachines - Vol. I," NASA CR 159635, 

July 1979. 

3. Watanabe, I., Arigh, I. and Fujie, K. , "A Study of the Flow 
Patterns in an Impeller Channel of a Radial Turbine," Trans. 
ASME, Journal of Engineering for Power, Oct. 1967, p. 467. 

4. Benson, R.S., Cartwright, W.G. and Hill, M. , "Analytical 
and Experimental Studies of Two-Dimensional Flows in a 
Radial Bladed Impeller," ASME Paper 71-GT-20, 1971. 

5. Katsanis, T. , "Fortran Program for Calculating Transonic 
Velocities on a Blade-to-Blade Stream Surface of a Turbo- 
machine," NASA TN D-5427, 1969. 

6. Wu, C.H. , "A General Theory of Three-Dimensional Flow in 
Sxdjsonic and Supersonic Turbomachines of Axial, Radial and 
Mixed Flow Types," National Advisory Committee for Aero- 
nautics, NACA TN 2604, 1952. 

7. Khalil, I. and Tabakoff, W. , "A Calculation Procedure for 
Viscous Flow in Turbomachines - Vol. II," NASA CR 159636, 
January 1980. 

8. Launder, B.E. and Spalding, D.B., "The Numerical Computations 
of Turbulent Flows," Computer Methods in Applied Mechanics 
and Engineering, Vol. 3, 1974, pp. 269-289. 

9. Wilcox, D.C. and Chambers, T.L. , "Critical Examination of 
Two-Equations Turbulence Closure Models," AIAA Paper No. 
76-352, 1976. 

10. Quemord, C. and Michel, R. , "Definition and Application of 
Means for Predicting Shear Turbulent Flows in Turbomachines," 
ASME Paper 76-GT-67, 1976. 

11. Dodge, P.R., "A Numerical Method for 2-D and 3-D Viscous 
Flow," AIAA Paper No. 76-425, 1976. 


21 


12. Roache, P.J./ Computational Fluid Dynamics ^ Hermosa Publishers, 
Albuquerque, New Mexico, 1972. 

13. Thompson, J.F., Thames, F.C. and Wayne, C., "Automatic 
Numerical Generation of Body-Fitted Curvilinear Coordinate 
System for Field Containing any Number of Arbitrary Two- 
Dimensional Bodies," Journal of Computational Physics, 

Vol. 15, 1974, pp. 299-319. 


22 





FIG. 1.1. HUB-SHROUD PROFILE WITH THE STkEAM SURFACE Sj. 
USED FOR BLADE-TO-BLADE ANALYSIS. 






PREDICTION CASE 2 


FIG. 1.3. NONDlfENSIONAL VELOCITY DISTRIBUTION AT 
DIFFERENT RADIAL LOCATIONS. 




HUB-SHROUD PROFILE WITH THE STREAM SURFACE, Sj, USED FOR 
BLADE-TO-BLADE ANALYSIS (COMPRESSOR CASE). 




PREDICTED 



EXPERIMENTAL 


FIG, 1.5. RELATIVE STREAMLINES FOR FLOW THROUGH 
RADIAL COMPRESSOR. 














INPUT FORM 















































Sadlal Blade 



HUB-SHROUD PROFILE WITH STREAM SURFACE Sj, USED FOR 
BLADE-TO-BLADE ANALYSIS. 
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IRID ORDERING SYSTEM IN THE TRANSFORM DOMAIN 
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